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A calculation for polar Kerr effect in high temperature cuprate superconductors 
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A mechanism is proposed for the tantalizing evidence of polar Kerr effect in a class of high temperature 
superconductors-the signs of the Kerr angle from two opposite faces of the same sample are identical and 
magnetic field training is non-existent. The mechanism does not break global time reversal symmetry, as in an 
antiferromagnet, and results in zero Faraday effect. It is best understood in a phenomenological model of bilayer 
cuprates, such as YBa 2 Cu 3 06 + 5 , in which intra-bilayer tunneling nucleates a chiral d-density wave such that 
the individual layers have opposite chirality. Although specific to the chiral d-density wave, the mechanism 
may be more general to any quasi-two-dimensional orbital antiferromagnet in which time reversal symmetry is 
broken in each plane, but not when averaged macroscopically. 


I. INTRODUCTION. 


The origin and nature of the pseudogap phase in the high- 
Tc c^ate superconductors still remains an unresolved prob¬ 
lem The pseudogap phase, which occurs in the under¬ 
doped regime of hole doping, and at temperature range T* > 
T > Tc, displays many interesting properties including vari¬ 
ous charge, spin, electron nematic, or current ordered states 
competing with superconductivit}!®!^. Recently a nonzero 
polar Kerr effect (PKE) has been observed in the pseudogap 
phase in a number of recent experimental^®!] but with un¬ 
usual characteristics. The effect measures the angle of rotation 
of linearly polarized light reflected from a medium at normal 
incidence and typically signals time reversal symmetry (TRS) 
breaking in the reflecting mediunP^. In a ferromagnetic ma¬ 
terial, the signs of the polar Kerr angle from two opposite sur¬ 
faces of the same sample are expected to be different. This 
is because the net magnetic moment points in the same direc¬ 
tion throughout the sample and hence if it points away from 
the sample on the top surface, it points into the sample on the 
bottom surface, see Fig .[2 Moreover, It should be possible to 
choose (or Train’) the direction of the net magnetic moment, 
and, in turn, the sign of the polar Kerr angle, by cooling the 
sample in the presence of a magnetic held. 

In contrast, in high-Tc superconductors it has been observed 
that the signs of PKE from the two opposite surfaces of the 
same sample are identical, and, moreover, the signal cannot be 
trained by magnetic held. To account for these puzzling ex¬ 
perimental observations, time-reversal invariant models with 
gyrotropic order were employed recentl}^^^^. However, the 
concept of gyrotropic order as an explanation for non-zero 
PKE in the cuprates were subsequently retractecP^®^ because 
it does not satisfy Onsager’s reciprocity principle in normal 
reflection that forbids a non-zero PKE in the absence of TRS 
breakin^^S®. 

Here we show that the observations in high-Tc can be un¬ 
derstood in the framework of a chiral d-density wave stat^ 
in the presence of interlayer tunneling, which is invariant un¬ 
der TRS in the bulk, but can still have a non zero PKE be¬ 
cause it is a property of the light reflected from the top surface 
which breaks TRS locally. Therefore PKE in this mechanism 
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FIG. I. Left panel: Ferromagnetic (FM) ordering where magnetic 
moments point in the same direction at each plane, thus resulting in 
the opposite signs of the Kerr angle from the top and bottom surfaces. 
Right panel: Orbital antiferromagnetic (OAF) ordering, where mag¬ 
netic moments switch sign on each plane, so the Kerr angle has the 
same sign from both surfaces. In contrast to FM, magnetic moments 
in OAF point out of the sample at the top and the bottom surfaces, 
provided the sample is cleaved between the bilayers. 


would be insensitive to the changing skin depth of the incom¬ 
ing light at the top surface, while PKE from a bulk order pa¬ 
rameter based description would yield a stronger effect for a 
longer skin depth. The chiral d-density wave state is deflned 
by adding a small dxy component to the dominant 
i.e. with the combined order parameter dxy + idx2-y2. The 
net order parameter breaks TRS at each CuO plane and results 
in a non-zero Hah conductivity The addition of a possi¬ 

ble dxy component could be a result of microscopic electronic 
interactions^U, or a structural transition that breaks the symme¬ 
try between the neighboring plaquettes. In the idx‘2-y2 state, 
by itself, spontaneous currents alternatingly circulate around 
plaquettes of the two-dimensional square lattice, thus preserv¬ 
ing the macroscopic TRS, but not any associated chirality. 

We establish that, for our present model, the angle of ro¬ 
tation due to one layer is cancelled by its neighbor, result¬ 
ing in zero Earaday rotation of the polarization plane of the 
transmitted light. However, since PKE is primarily a sur¬ 
face phenomenon, where the light reflected from the top (or 
bottom) surface at normal incidence changes its plane of po¬ 
larization, there can be a non-zero PKE. Eurthermore, since 
bilayer cuprates usually cleave through the reservoir layers 
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FIG. 2. (Color online) Left: Band dispersion as a function of k = 
(/cxjTt) for the two states in a bilayer: d -\- idjd -\- id (blue) and 
d + id/d — id (red). We utilized t' = 0.33t, Wo = .33t, Ao = 
O.OlWo/2, and t± = 0.5t all energies measured in the unit of t. 
Right: The ground state energy versus hole doping (from 0.08 to 
0.18) indicating that ihQ d + id/d — id state (in red) has lower energy 
in a bilayer for any given value of hole doping. 


separating the CuO bilayers, the magnetizations at the top 
and bottom surfaces should point opposite to each other (see 
Fig. [T]), giving rise to the same sign of the Kerr angle. Fi¬ 
nally, since the system as a whole is an OAF, coupling to a 
small external magnetic field should be small, resulting, most 
likely, in a small or non-existent magnetic field ‘training’ ef¬ 
fect. Importantly, to the best of our knowledge, the scenario 
presented here is the only one consistent with all the puz¬ 
zling phenomenology seen in the recent PKE experiments in 
cuprates. 

It has been argued that much of the phenomenol ogy of the 
cuprates in the underdoped regime can be unifiecpEIlIllI by 
making a single assumption that the ordered -density 

wave (DDW) state is responsible for the pseudogap. More¬ 
over, an extensive Hartree-Fock calculation for idx2_y2 state 
has recently been carried ou® So far, evidence of mag¬ 
netism arising from (i-density wave in neutron or NMR mea¬ 
surements has been controversial. However, the success of 
the present phenomenological model in explaining PKE must 
speak in favor of the suggested order parameter. 

This paper is divided as follows: in Section II, we introduce 
the chiral (i-density wave state order parameter, and calculate 
the anomalous Hall conductivity of a single layer. In Sec¬ 
tion III, we discuss the problem of light propagation through 
a single cuprate layer, and then calculate the Kerr and Fara¬ 
day responses through the bilayer system in Section IV. We 
conclude in Section V. 


II. CHIRAL DDW WITH INTERLAYER TUNNELING. 

Consider a combination of the density waves, d^y + 
idx2_y2 , such that the net order parameter is 

(^k+Q,a^k,/3) ^ ['^^k “ (1) 

where cj^ ^ is the electron creation operator of momentum k 
and spin a, and Q = (tt, tt) is the density wave vector. W\^ 


and Ak correspond to idx2_y2 and dxy respectively, defined 
as FFk = ^{coskx — cosky), Ak = Aq sin kx sin ky. We 
consider a bilayer system where the idx2-y2 component of the 
order parameter, i.e., W\^ may or may not switch sign between 
the two layers. The four component mean field Hamiltonian 
for the system in the basis = (Ck^, Ck+Q, Ck^, ^^k+q) t^kes 
the following form 


H{k) = 


/ ^k ^ik ^±k 0 \ 

^Ik ^k+Q 0 t±k+Q I 

^±k 0 Ck ^2k I 

V 0 t±k+Q 921 ^ Ck+Q / 


( 2 ) 


where Ck is the energy dispersion for a two-dimensional 
square lattice. 


Ck = —2f(cos kx + cos ky) + cos kx cos ky^ (3) 

where t and t' are the nearest and next-nearest hopping inte¬ 
grals in the tight-binding Hamiltonian, 5^1 k = ^kUk — A\^ and 
^±k = t_\_{coskx — cos ky )‘^/4 describes the tunneling be¬ 
tween the two layer® appropriate for tetragonal systems. The 
superscript ( 1 , 2 ) on the electron operator in is the layer in¬ 
dex. Note that ^ik = ^2k represents 3. d id/d id bilayer 
configuration and = ^2k is Sid^ id/d —id configuration. 
We find that when = ^2k the system is energetically more 

favorable than the case when ^ik = ^2k- This is observed 
by diagonalizing the Hamiltonian and obtaining the ground 
state energy for a given doping concentration, as displayed in 
Fig. The d ^ id state spontaneously breaks time-reversal 
symmetry (T) as well as the in-plane refiection symmetry 
about the principal axes and exhibits anomalous Hall effect 
with a non-zero value of the Hall conductivity cFxy. However, 
the value of axy reverses sign for the d — id state. Thus the 
ground state of the bilayer breaks TRS in each plane, but since 
the inversion symmetry (P) about the mid point between the 
planes is also broken, the product PT is conserved, allowing 
the system to have a nonzero polar Kei r effe ct despite conserv¬ 
ing the global TRS and being an OAF*^^^. The magnetoelec¬ 
tric effect and PKE in a nothe r antiferromagnet Cr203 were 
predicted theoretically in I 3 S 1 , and subsequently observed in 
experiment®®. 

The Hall conductance of a single layer described by a 
d P id mean field Hamiltonian can be calculated u sing the 
formalism of linear response theory and Kubo formul®. The 
two-component mean field Hamiltonian describing di d -f id 
density-wave state in the = (cj^, basis is given by: 




ek g\^ ^ 

^k+Q ) 


( 4 ) 


At a finite frequency cc and in the limit q ^ 0 , the anomalous 
Hall conductivity at any finite temperature is given b}® 


^ 2 ^ r de p(k)/(p+(k))-/(p_(k)) 
h J { 2 n)^ w{]P)[z — 2 w{]P)][z P 2 w{k.)]^ 


where P(k) = 4 :tAoWo{sin^ ky + cos‘^ ky sin^ kx) is the 
Berry curvature, u;(k) is the modulus of a three compo¬ 
nent vector w(k) = [-Ak, -VFk, (^k - ek+Q)/2], / is the 
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Fermi distribution function, /i is the chemical potential, and 
2 ; = cc + i(5, with S a positive infinitesimal. F^±(k) = 
(^k + ek+Q)/2 ± u;(k) — /i describe the two energy bands 
obtained by diagonalizing the Hamiltonian in Eq. Q. The 
sign of is determined by the sign of the product AoVEo. 
so the d±id states have opposite signs of axy. 



\/\/\/ 

\\\ 


z=-h ^ 


III. TRANSMISSION AND REFLECTION OF LIGHT 
FROM A SINGLE LAYER. 


We now study propagation of an electromagnetic wave 
through a layered system with chir al DDW using standard 
electrodynamics formalisnP^I^^I^^I^. First we consider an 
electromagnetic wave incident normally on a single two- 
dimensional layer of a material in the xy plane. The electric 
field components of the wave in a medium are given by 


E = 


K 

El 


-h 


El 

EL 


( 6 ) 


El and EL are the transmitted components of right and left 
circularly polarized (CP) light respectively and similarly El 
and EL are the refiected components, and k is the wavevec- 
tor. The corresponding magnetic field components can be 
found using Maxwell’s equation, k x E = ujH. The com¬ 
ponents of the electromagnetic field satisfy standard electro¬ 
dynamic boundary conditions at the material layer, which we 
assume is located at z = h: E^h = {Hyh — H<h)y = 
—4:7r{dE)x, {H^h — H<h)x = ^T^{dE)y. Note that a is a ten¬ 
sor, but ^ is a two component vector. We consider only the 
Hall components of the surface current and neglect the pres¬ 
ence of (7xx\ however, we have confirmed that the presence 
of a finite cr^x due to the chiral DDW Fermi pockets does not 
change our results qualitatively. We define the scattering ma¬ 
trix which relates the incoming and the outgoing electric field 
components, O = SI, where the outgoing wave O and the 
incoming wave I are given by. 



” ITT 

^+,<h 


\E<h^ 

0 = 

E>h 

and I = 

rpr 


[E>h\ 




and the scattering matrix S is. 




r R++ i?+_ t;+ t;_ 1 

n V ■ 


1 

+ 

1 

1 

+ 

1 

r w 


T++ T+_ R\_ 



. r_+ r__ R'__ . 


This scattering matrix S describes reflection and transmission 
of electric field components from the top surface of the slab. 
We have also defined in Eq. (§ two-component matrices IZ, 
TL T and IZL whose components are given by the corre¬ 
sponding block entries. Matching the boundary conditions at 


FIG. 3. (Color online) Schematic diagram showing multiple reflec¬ 
tions and transmissions through the top and bottom layers, which we 
use to calculate the Kerr and Faraday angles. Note that we have as¬ 
sumed normal incidence for the incoming light in our calculations. 


z = h,wQ find that 

^ (1 - - (dTTCr^^)^ + iSTTCTxy) 

(1 + nf + {Airaxy)^ 

^ (1 - - {Airaxy)^ - iSiraxy) 

(1 -h nf + (47rcr^^)2 

edk>-k<)h ^2 (^1 ^ ^ 

(1 + nf + {4:7r(Jxy)‘^ 

^ (2 (1 + n^) — iSTTCTa^y) 

(1 + nf + {^'KCFxyY 

where n is the refractive index of the medium. R' _, 

T^+,and TL_ can be obtained in a similar fashion. The other 
components of S which couple right and left CP components 

i.e. _, TL_^ and so on, all vanish. We note that when a^y 7 ^ 

0, 7 ^ R _and T++ 7 ^ T_which is a signature of 

broken time-reversal symmetry. 


IV. POLAR KERR AND FARADAY EFFECTS IN BILAYER 
CHIRAL DDW. 

To discuss scattering from the bilayer, we consider two such 
interfaces ai z = Eh and z = —h as depicted in Fig.[^ Since 
even in the presence of the interlayer coupling t±{k) the sys¬ 
tem brea ks P and T while PT is conserved, allowing a non¬ 
zero PKri®3 in the following we will ignore t_\_{k) for sim¬ 
plicity, expecting it to modify our results only quantitatively. 
The Hall conductivity a^y reverses it sign at the bottom layer 
at z = —h. One can then appropriately define the scatter¬ 
ing matrix elements for the bottom layer taking into account 
the opposite sign of the Hall conductivity and the position of 
the bottom plane to be —h instead of h. We denote the two- 
component matrices defined in Eq. ^ for the top layer by the 
subscript T and by subscript B for the bottom layer. Thus 
refiection and transmission through the bilayer as whole are 
described by tensors R and T given by, 

R = 7Zt + R'tRbL^ ~ R!tIZb) ^7t 
T = TB{^-n'TnBETT (9) 

We now switch basis from CP light to linearly polarized (LP) 
light for convenience of the following discussion. Denoting 
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the electric field of the light incident on the sample by Ej, 
Er = HEi and Er = TEi give the refiected and the trans¬ 
mitted electric fields. When linearly polarized light is incident 
on the sample, the Kerr and Faraday angles are determined by 
the difference between right and left CP light: 

Of = ^(arg[£’+] - arg[£’“]) 

Ok = ^{arg[E+] - aig[E^]), (10) 


where E^ ^ for Er^r = [Er^r, E'^ ^]. For 

the bilayer system discussed above, the R has non-zero off 
diagonal elements (in LP basis) and T is diagonal, which is a 
clear signature of a non-zero Kerr response and the absence of 
the Faraday effect. (We do not state the analytic expressions 
for these matrices here, as they are too cumbersome.) 

We now make a rough estimate for the polar Kerr angle 
for a bilayer system using Eqs. (|^ and (1^). Measuring all 
the energies in units of t we use t' = 0.3, jj. = —.9, n ^ 
the interlayer distance 2= 3.2 A, the strength of the 
idx 2 _y 2 component of the order parameter Wo{p) = 0.1(1 — 
p/Pc), where p is the hole doping concentration and Pc = 
0.17, Ao = O.OllFo/2 and the frequency of measurement 
uj = 1500 nm. In Fig.|^ we have plotted the polar Kerr angle 
Ok as a function of hole doping and we obtain a non-zero Ken- 
angle of the order of 100 nrad. The estimated Faraday angle 
from our formalism turns out to be zero, again from Eqs. 0 
and ( p^ . Since the chiral DDW with interlayer tunneling is an 
OAE, the angle of rotation of the plane of polarization of light 
due to one layer is cancelled by its neighbor, resulting in zero 
Earaday rotation of the transmitted light. However, since PKE 
is primarily a surface phenomenon, where the light refiected 
from the top surface changes its plane of polarization, there 
is a non-zero PKE. Eurther, since the magnetizations at the 
top and bottom surfaces should point opposite to each other 
(see Eig. [T]), the two surfaces give rise to the same sign of 
the Kerr angle. Einally, since the system as a whole is an 
OAE, coupling to a small external magnetic field should be 
small, leading to small or non-existent ‘training’ effect. It is 
important to note that all of these conclusions are consistent 
with the phenomenology of the recent PKE measurements in 
the cuprates. 


V. CONCLUSIONS. 

To conclude we considered the chiral DDW state in a bi¬ 
layer where the sign of the idx 2 _y 2 component of the order 
parameter changes between the layers which is an energeti¬ 


cally more favorable configuration. This also leads to the re¬ 
versal of sign of axy in the bottom layer, thus breaking inver¬ 
sion symmetry. The calculations presented here are consis¬ 
tent with the unusual PKE observed in high-Tc materials. Our 
calculations, although applied here specifically to the chiral 
DDW state, are more generally valid for any OAE with TRS 
broken at each plane. In Ref.® similar ideas were applied to 
a tilted loop current model. In addition, the ideas presented 
here also apply to the bi-axial density wave recently seen in 
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FIG. 4. Estimated Kerr angle in nrad as function of hole doping 
p. The strength of the id^ 2 _y 2 is assumed to vary with doping as 
W{p) = 0.1(1 — p/pc) eV, where pc is chosen to be 0.17. The am¬ 
plitude Ao of the dxy component is assumed to be 1% of Wo. While 
a non-zero PKE is a robust consequence of our model the precise 
values of Wo, Ao, Ok in this figure are for illustrative purposes only. 



the pseudogap phasd^^Jllj [f acconipanied by sponta¬ 

neous currents®E2 However, the theories®E3 are currently 
formulated for a single layer, and it remains to be seen whether 
they can be generalized to a multilayer model with alternating 
sign of axy similar to the present work. 

Another important class of high-Tc materials is single layer 
compounds, such as Bi-2201 (Bi2+xSr2-xCu08+(5) and Hg- 
1201 (HgBa 2 Cu 04 +< 5 ). Although the detailed results are not 
yet published, it is known that such materials also show simi¬ 
lar PKE, as discussed herJ^. At the level of order parameter 
symmetry, there is no difference, in the sense that one can eas¬ 
ily envision CuO-layers alternating between d-\-id and d — id. 
In addition, recent X-ray measurements indicate that the unit 
cell in the c-direction is doubled, bringing it closer to the bi¬ 
layer problem. Until PKE measurements in single layer ma¬ 
terials are published in detail, it is probably prudent to refrain 
from further speculations. 

Acknowledgement. G.S. and S.T. are supported by AEOSR 
(EA9550-13-1-0045). P.G. is supported by JQI-NSE-PEC. 
S.C. was supported by a grant from NSE-DMR-1004520. 


^ M. R. Norman, D. Pines, C. Kallin, Adv. Phys. 54, 715 (2005). 

^ B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida, J. Zaanen, 
Nature 518, 179 (2015). 

^ C. M. Varma, Phys. Rev. B 55, 14554 (1997); ibid 73, 155113 
(2006). 


S. Chakravarty, R. B. Laughlin, D. K. Morr, C. Nayak, Phys. Rev. 
B 63, 094503 (2001). 

^ S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, J. M. 
Tranquada, A. Kapitulnik, and C. Howald, Rev. Mod. Phys. 75, 
1201 (2003). 








5 


^ S. Chakravarty, Phys. Rev. B 89, 087101 (2014). 

^ T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. Berthier, W. N. 
Hardy, R. Liang, D. A. Bonn, and M.-H. Julien, Nature 477, 191 
( 2011 ). 

^ T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. Berthier, P. L. 
Kuhns, A. P. Reyes, R. Liang, W. N. Hardy, D. A. Bonn, and M.- 

H. Julien, Nature Communications 4, 2113 (2013). 

^ N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, J.-B. 
Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, L. Taillefer, 
Nature 447, 565 (2007). 

S. E. Sebastian, N. Harrison, E. Palm, T. P. Murphy, C. H. Mielke, 
R. Liang, D. A. Bonn, W. N. Hardy, G. G. Lonzarich, Nature 454, 
200 (2008). 

J. Chang, E. Blackburn, A. T. Holmes, N. B. Christensen, J. 
Larsen, J. Mesot, R. Liang, D. A. Bonn, W. N. Hardy, A. Wa- 
tenphul, M. v. Zimmermann, E. M. Eorgan, S. M. Hayden, Nature 
Phys. 8, 871 (2012). 

G. Ghiringhelli, M. Le Tacon, M. Minola, S. Blanco-Canosa, 
C. Mazzoli, N. B. Brookes, G. M. De Luca, A. Erano, D. G. 
Hawthorn, E. He, T. Loew, M. Moretti Sala, D. C. Peets, M. Sal- 
luzzo, E. Schierle, R. Sutarto, G. A. Sawatzky, E. Weschke, B. 
Keimer, L. Braicovich, Science 337, 821 (2012). 

J. Xia, E. Schemm, G. Deutscher, S. A. Kivelson, D. A. Bonn, W. 
N. Hardy, R. Liang, W. Siemons, G. Koster, M. M. Eejer, and A. 
Kapitulnik, Phys. Rev. Lett. 100, 127002 (2008). 

A. Kapitulnik, J. Xia, E. Schemm, and A. Palevski, New J. Phys. 
11, 055060 (2009). 

R.-H. He, M. Hashimoto, H. Karapetyan, J. D. Koralek, J. P. Hin¬ 
ton, J. P. Testaud, V. Nathan, Y. Yoshida, H. Yao, K. Tanaka, W. 
Meevasana, R. G. Moore, D. H. Lu, S.-K. Mo, M. Ishikado, H. 
Eisaki, Z. Hussain, T. P. Devereaux, S. A. Kivelson, J. Orenstein, 

A. Kapitulnik, and Z.-X. Shen, Science 331, 1579 (2011). 

H. Karapetyan, M. Hiicker, G. D. Gu, J. M. Tranquada, M. M. 
Eejer, J. Xia, and A. Kapitulnik, Phys. Rev. Lett. 109, 147001 
( 2012 ). 

H. Karapetyan, J. Xia, M. Hiicker, G. D. Gu, J. M. Tranquada, M. 
M. Eejer, and A. Kapitulnik, Phys. Rev. Lett. 112, 047003 (2014). 

B. 1. Halperin in The Physics and Chemistry of Oxide Super¬ 
conductors, Vol. 60, edited by Y. lye and H. Yasuoka (Springer, 
Berlin, Heidelberg, 1992). 

P. Hosur, A. Kapitulnik, S. A. Kivelson, J. Orenstein, and S. 
Raghu, Phys. Rev. B 87, 115116 (2013). 

J. Orenstein and J. Moore, Phys. Rev. B 87, 165110 (2013); see 
also the comment by S. Chakravarty, Phys. Rev. B 89, 087101 
(2014). 

V. P. Mineev, Phys. Rev. B 88, 134514 (2013). 

22 S. S. Pershoguba, K. Kechedzhi, and V. M. Yakovenko, Phys. Rev. 
Lett. Ill, 047005 (2013). 


V. P. Mineev and Yu. Yoshioka, Phys. Rev. B 89, 139902(E) 
(2014). 

S. S. Pershoguba, K. Kechedzhi, and V. M. Yakovenko, Phys. Rev. 
Lett. 113, 129901(E) (2014). 

P. Hosur, A. Kapitulnik, S. A. Kivelson, J. Orenstein, S. Raghu, 

W. Cho, and A. Pried, Phys. Rev. B 91, 039908(E) (2015). 

N. P. Armitage, Phys. Rev. B 90, 035135 (2014). 

A. D. Pried, Phys. Rev. B 90, 121112 (2014). 

A. Kapitulnik, Physica B 460, 151 (2015). 

W. Cho, S. A. Kivelson, arXiv: 1507.03555 (2015). 

S. Tewari, C. Zhang, V. M. Yakovenko, and S. Das Sarma, Phys. 
Rev. Lett. 100, 217004 (2008). 

C. Nayak, Phys. Rev. B 62, 4800 (2000). 

S. Chakravarty et al, Int. J. Mod. Phys. B15, 2901 (2001); A. 
Kopp and S. Chakravarty, Proc. SPIE, Vol. 5932, 593219 (2005) 
I arXiv:cond-mat/0507574j. 

S. Tewari et al., Phys. Rev. B 70, 014514 (2004). 

S. Chakravarty et al.. Nature 428, 53 (2004). 

R. B. Laughlin, Phys. Rev. Lett. 112, 017004 (2014); Phys. Rev. 
B 89, 035134 (2014). 

1. E. Dzyaloshinskii, Phys. Lett. A 155, 62 (1991); 1. Dzyaloshin- 
skii and E. V. Papamichail, Phys. Rev. Lett. 75, 3004 (1995). 

G. S. Canright and A. G. Rojo, Phys. Rev. Lett. 68, 1601 (1992); 
Phys. Rev. B 46 14078 (1992). 

I. E. Dzyaloshinskii, Zh. Eksp. Teor. Piz. 37, 331 (1959) 

R. M. Hornreich and S. Shtrikman Phys. Rev. 171, 1065 (1968) 

D. N. Astrov, Zh. Eksp. Teor. Piz. 38, 984 (1960) 

V. J. Polen, G. T. Rado, and E. W. Stalder, Phys. Rev. Lett. 6, 607 
(1961) 

G. T. Rado and V. J. Polen, Phys. Rev. Lett. 7, 310 (1961) 

B. B. Krichevtsov, V. V. Pavlov, R. V. Pisarev, and V. N. Gridnev, 

J. Phys. Cond. Matt. 5, 8233 (1993) 

B. B. Krichevtsov, V. V. Pavlov, R. V. Pisarev, and V. N. Gridnev, 
Phys. Rev. Lett. 76, 4628 (1996) 

W. K. Tse and A. H. MacDonald, Phys. Rev. Lett. 105, 057401 

( 2010 ). 

W. K. Tse and A. H. MacDonald, Phys. Rev. B 84, 205327 (2011). 

H. Kezuka, et al., Physica C: Superconductivity 185, 999 (1991). 
V. M. Yakovenko, Physica B 460, 159 (2015). 

Y. Wang and A. V. Chubukov, Phys. Rev. B 90, 035149 (2014). 

D. P. Agterberg, D. S. Melchert, and M. K. Kashyap, Phys. Rev. 
B 91, 054502 (2015). 

Y. Wang, A. V. Chubukov, and R. Nandkishore, Phys. Rev. B 90, 
205130 (2014). 

M. Gradhand, I. Eremin, and J. Knolle, Phys. Rev. B 91, 
060512(R) (2015). 

A. Kapitulnik, private communication. 


23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 


